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The scattering processes of exotic atoms in excited states from hydrogen such as 
elastic scattering, Stark transitions and Coulomb de-excitation are studied within a 
close coupling approach. The vacuum polarization and the strong interaction shifts 
of ns-states (in case of hadronic atoms) are taken into account. The differential and 
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integral cross sections of the above processes are calculated to use them as the input 
in cascade calculations. The effect of closed channels on the scattering processes is 
investigated. 
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1 Introduction 

The exotic hydrogen-like atoms are formed in highly excited states, when negative particles 
7r~, K~ ...) are stopped in hydrogen. The deexcitation of exotic atoms proceeds via 
many intermediate states until the ground state is reached or a nuclear reaction takes place. 
Despite a long history of theoretical and experimental studies the kinetics of this atomic 
cascade is not yet fully understood. The present experiments with the exotic hydrogen-like 
atoms addresses a number of fundamental problems using precision spectroscopy methods, 
the success of which relies crucially on a better knowledge of the atomic cascade. The 
experimental data are mainly appropriate to the processes at the last stage of the atomic 
cascade (X-ray yields and the products of the weak or strong interaction of the exotic 
particle in the low angular momentum states with hydrogen isotopes). So the reliable 
theoretical backgrounds on the processes both in low-lying and in highly excited states are 
required for the detailed and proper analysis of the experimental data. 

In this paper we present the ab initio quantum-mechanical treatment of non-reactive 
scattering processes of the excited exotic hydrogen atom in collisions with the hydrogenic 
atom in the ground state: 

{aX) n! i + {be~) ls -> (aX)tfjt + (be~) ls , (1) 

elastic scattering^' = n, I' — I), Stark transitions (n' = n, I'^l), and Coulomb deexcitation 
(CD) in' <n). Here (a, b) = (p,d,t) are hydrogen isotopes and X = , ir~ , K~ , p; (n,l) are 
the principal and orbital quantum numbers of exotic atom. While the deexcitation processes 
are obviously essential for the atomic cascade, the role of the collisional processes preserving 
the principal quantum number n is also very important. The Stark transitions affect the 
population of the nl sublevels and together with the elastic scattering they decelerate the 
exotic atoms thus influencing their energy distribution during the cascade. 
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Starting from the classical paper by Leon and Bethe pQ, Stark transitions has been 
treated in the semiclassical straight-line-trajectory approximation (see [2] and references 
therein) . 

The fully quantum-mechanical treatment of the elastic scattering and Stark transitions 
based on the adiabatic description was given for the first time in [31 0]. Recently [5], the 
processes have also been studied in a close-coupling approach with using dipole approxima- 
tion for interaction potential and taking the electron screening effect into account by the 
model. The cross sections calculated in approaches [SI HIE] are in good agreement. 

Concerning the CD process, the situation is much less defined, especially for low n. 
The first work on the CD process was performed by Bracci and Fiorentini [6] in frame 
of the semiclassical approach with some additional approximations. In the following nu- 
merous papers (see [7J [8] and references therein) the CD process is considered within the 
asymptotic approaches using the adiabatic hidden crossing theory. The CD calculations 
were also performed in the classical-trajectory Monte-Carlo (CTMC) approach^. While 
the Coulomb deexcitation cross sections obtained in CTMC approach are in fair agreement 
with the semiclassical ones of Bracci and Fiorentini [6] , the more elaborated advanced adi- 
abatic approach (AAA) [TJ [8] gives too small CD cross sections to explain the experimental 
data JTU]. The reasons of such a strong discrepancy are not clear. One can only assume 
that the semiclassical model of Bracci and Fiorentini as well as the CTMC approach are 
not valid for low-lying states and at low energies. 

The processes (1) have been treated recently by authors in a unified manner in the 
framework of the close-coupling (CC) approach (see for detail [HI [12j [131 EE])- The dif- 
ferential and integral cross sections for the processes (1) have been calculated for muonic, 
pionic and antiprotonic hydrogen atoms in excited states with n = 2 — 14 and in a kinetic 
energy range relevant for cascade calculations. The energy shifts of the ns states due to 
vacuum polarization and strong interaction (for hadronic atoms) are included in the close- 
coupling method. This approach allows to obtain the self-consistent description of all the 
processes (1) and is free from the additional approximations used in previous studies. The 
calculated differential and integral cross sections presented in this paper mainly refer to the 
CD process to illustrate some of our new results obtained quite recently. 

All open channels corresponding to the exotic atom states with n <no (no is principal 
quantum number in the entrance channel) have been included in the CC calculations. The 
effect of the closed channels with n > n was also studied and will be discussed below. 
(Throughout the whole paper the cross sections are given in atomic units.) 

2 Close-coupling approach 

The Hamiltonian of the system ((X~p) n i + Hi s ) (after separation of the cm. motion) is 
given by 



where m is the reduced mass of the system, R is the radius vector between the cm. of the 
colliding atoms, p and r are their inner coordinates. The interaction potential, V(r,p, R), 
is a sum of four Coulomb pair interactions between the projectile atom and the target 
atom particles. h ex (p) and hn(r) are the hydrogen-like Hamiltonians of the free exotic and 
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hydrogen atom, whose eigenfunctions together with the angular wave function 1la(R) of 
the relative motion form the basis states, |ls, nl, L : JM), with the conserving total angular 
momentum (JM) and parity tt = (— l) l+L . In the present consideration we use the "frozen" 
electron approximation. The CC approach can be extended in a straightforward manner to 
include the target electron excitations. The total wave function of the system are expanded 
in terms of the basis states as follows 

* J E M *(r, p, R) = R- 1 J2 G J J L (R)\U, nl, L : JM). (3) 

nlL 

The expansion (3) leads to the close-coupling second order differential equations for the 
radial functions of the relative motion, G^J L (R), 

(^ + k n- L{L m l) ) G ^( R ) = Z™J2 W^ VL , ML (R) G J J VL ,(R). (4) 

^ ' n'l'L' 

The channel wave number is defined as k 2 = 2m(E cm + E no i — E n i), where E cm and (n l ) 
are the energy of the relative motion and the exotic atom quantum numbers in the entrance 
channel, respectively. The bound energy of the exotic atom, E ni = e n \ + Ae n i, includes the 
eigenvalue of h ex (p), e n i, and the energy shift, Ae n i = Ae^ + Ae^f, due to the vacuum 
polarization and strong interaction (in case of adronic atom). Hereafter, the energy E cm 
will be referred to £„^o i n the entrance channel (we assume here that Ae n ^ = 0) 
The matrix elements of the interaction potential, V(r, p, R), 

Wn'vv,niL{ R ) = n ' 1 '» L ' : JM\V\U, nl, L : JM) (5) 

are obtained by averaging it over the electron wave function of the ls-state and then 
applying the multipole expansion. The integration over (p, R) reduces the matrix elements 
(5) to the multiple finite sum. 

At fixed E cm the coupled differential equations (4) for the given J and tc values are solved 
numerically by the Numerov method with the standing-wave boundary conditions involving 
the real symmetrical .fT-matrix related to T-matrix by the equation T = 2iK/(I — iK). 
All open channels corresponding to the exotic atom states with n < n have been included 
in the close-coupling calculations. The effect of closed channels with n > n was also 
investigated and will be discussed below. In the next sections we present and discuss the 
following total cross sections of the scattering processes: the partial cross section a J 
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the total cross section of the nl — > n'l' transition o n i~>n'v = cr nZ^n'/'(-^'); an d the l- 
averaged cross section a nn i = 1/n 2 ^ i ,(2Z / +l)cr n ^ n /;/(i?), and also the analogous differential 
cross sections. 



3 The elastic scattering and Stark transitions 

The total cross sections calculated in adiabatic [1] and present CC approaches 

whole in good agreement at energies E cm > 1 eV . The angular distributions obtained 
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in the present CC and adiabatic [3] calculations coincide in the region of the diffraction 
maximums but demonstrate significant differences in backward hemisphere (see the left side 
of Fig.l). The typical angular distributions of the Stark nl — » nl' transitions for n = 5 are 




Figure 1: The /-averaged differential (elastic and Stark) cross sections (left): adiabatic 
model [3 J (dashed lines) and present CC (solid lines). The /-averaged differential Stark cross 
sections for (//p)s; + H collisions vs. cms scattering angle 9 cm at E cm = 1 eV (right). 

shown in Fig. 1 (right). It is well known [SJ [5] that cross sections of these processes are 
similar to the diffraction scattering (at energies more 1 eV) with a strong forward peak 
enhancing with increasing energy and a set of maxima and minima. While the elastic cross 
sections always have a strong peak at 6 cm = 0, the first maximum position in the Stark 
cross section depends on the Al = \l — l'\ value. In particular, for Al = 1 this maximum is 
at finite scattering angles as it is also remarked in [5]. 




E ,eV 



Figure 2: The /-averaged Stark cross sec- 
tions for (pp) n + H collisions (present - 
solid line, semiclassical model [5] for n = 
8 - triangles) 




Figure 3: The cross sections <J n i^,nV for 
collisions (pp) n= s + H. The dashed and 
solid lines connect the points correspond- 
ing to the calculations both with and 
without the ns-state energy shifts, re- 
spectively. 



Some of our CC results for (pp) n + H collision are shown in Figs. 2,3: the E-dependence 
of /-averaged Stark cross sections for n = 2 — 14 (Fig. 2) and cross sections of the nl — > nl' 
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transitions for n = 8 at E cm = 2 eV (Fig. 3). It is seen from Fig. 3 that ns-state energy 
shift due to strong interaction leads to essential suppression of both Stark ns — > nl > 
and elastic np — > np transitions (we used the value e ls = 721 eV for the ls-state shift and 
£ns — tis/n 3 for ns-states) at energies compared with the shift value. The same effect is 
observed for all hadronic atoms. 

4 Coulomb de-excitation 

The nature of the CD process is quite different from the one of the elastic or Stark processes. 
In contrast to elastic (Stark) scattering CD process is accompanied by the large energy 
release (tens and hundreds eV) and occur at smaller distances, so the details of the short- 
range interaction are more important for the treatment of CD process than for elastic 
processes. This difference between elastic and deeply inelastic processes is illustrated by 



Figure 4: The partial cross sec- 
tions (in a.u.) of elastic scat- 
tering a^ 7 and Coulomb de- 
excitation a 76 for (np) n =7 + H 
collisions versus the total angu- 
lar momentum J at the ener- 
gies: 0.1 eV (dashed), 2 eV (dot- 
ted) and 50 eV (solid). 
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Fig. 4 where the J dependence of the partial-wave 1-averaged cross sections <j 3 nn , for n = 7 
at three fixed energies 0.1, 2, and 50 eV is shown. It is seen that a substantial part of the 
CD cross section (cr/ 6 ) comes from the partial waves with rather a low J in contrast to the 
elastic (Stark) process. 

In spite of the value of the total CD cross section constitutes about few per cent of the 
total elastic cross section it is incorrect to treat CD in the framework of the perturbation 
theory. In each significant for CD partial wave the value of the CD cross section is com- 
parable with the elastic cross section. The CD process is determined by the short-range 
behaviour of the wave function which changes when new channels are included in calcula- 
tion. Therefore, to calculate the transition n — > n — 1 in a proper way it is impossible to be 
restricted with the two- level approximation (n and n — 1) and the states with other nearest 
n should be involved. 

We studied the dependence of the results on the number of included channels and found 
that the inclusion of the channels with n < n — 1 leads to a strong suppression of the main 
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n — > n — 1 transitions in comparison with the two-level CC approximation and due to this 
the total CD cross section is also suppressed [TTJ. So in all our CD calculations we included 
all the open channel with n < uq. It should be noted that all the previous calculations of 
CD realized within semiclassical or adiabatic approaches used a two-level approximation. 
It is obvious that two-level approximation is not absolutely suitable for the treatment of 
transitions with An > 1. In contrast to the elastic scattering and Stark transitions where 
the "dipole" approximation (t m = 1) and even more rough dipole potential (used in [5]) 
gives reasonable results (at not too low energies), in case of the deeply inelastic process 
such as CD the full interaction must be used as it is clear from the present study. 

4.1 Muonic atoms 

It is commonly believed [SI El E] that the CD cross sections at low energies behave like 
1/E. In order to reveal more explicitly the distinction from the 1/E behaviour the present 
Z-averaged CD cross sections multiplied by energy are shown in Fig. 5 in comparison with 
the results of the SC model [6] for n = 3, 5, 7, 9 and CTMC calculations [9] for n = 9. As 



Figure 5: The cross sections of 
Coulomb de-excitation (multi- 
plied by E cm ) for {jip) n + H col- 
lisions calculated in the present 
CC method (solid lines) in com- 
parison with the SC [6j (dashed) 
and CTMC [9] results (dotted). 
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it is seen, the energy dependence of the CC cross sections in the region 5 > E > 1 eV, as 
a whole, is in a qualitative agreement with the results [5| and [9J. At lower energies the 
present CC results reveal ~ 1/ \fE dependence in accordance with the Wigner threshold 
law (the similar behaviour is seen in the CTMC results [9] for n — 9) and in disagreement 
with 1/E dependence obtained in the SC model [6] and AA approach [3 [8]. 

The distribution over the final states n' is strongly different from the SC results [6] as 
illustrated in Fig. 6. The CC calculations predict that the transitions with An > 1 are 
strongly enhanced as compared with the results of the two- level approaches [H [7J [8] . The 
An > 1 transitions make up a substantial fraction (16% - 37%) of the total CD cross section 
for n > 4. 

4.2 Differential cross sections 

The angular distributions of CD was calculated for the first time in [13] . Earlier in the cas- 
cade calculations the angular distributions of the CD process are presumed to be isotropic. 




6 



0.1: 



T v-" 

c 

3 

o.ow— 



O.OOljjp" 



-0 



T 



♦ 

.0 



1)1 













♦ v 






,.o 






n 








▼ ▼ 0=5, SC 






V-V n = 5 ' cc 




A-A n=6, CC 




• 


♦»•♦ n=7, SC 
O-O n=7, CC 
□-□ n=8, CC 
• • 0=9, SC 
G-O n=9, CC 



0.1 



0.01 



-- 0.001 



Figure 6: Dependence of the 
CD cross sections on the final 
principal quantum number n' 
for the different initial n in the 
(fip)n + H collisions at E = 1 
eV. The dashed and dotted lines 
connect the points calculated in 
the present CC and in the semi- 
classical (SC) [6J approaches, re- 
spectively. 



The calculated cross sections for individual nl — > n'V transitions with An = 1 and 2 at 
energy E cm = 1 eV are shown in Fig. 7. In Fig. 8 the /-averaged cross sections for the 6 — > 5 
transition at different values of the relative energy from 0.01 up to 15 eV are presented. We 
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Figure 7: Differential CD cross sections for the individual transitions with An = 1 (left) 
and with An = 2 (right) for n = 5 at E cm = 1 eV. 

found that the angular distributions both of the individual and /-averaged cross sections 
(excluding very low energies) are far from isotropic: as a whole the scattering at 9 cm < 60° 
and 8 cm > 120° is noticeably enhanced. The cross sections for ns —> n's transitions (see 
Fig. 7) have (as for elastic scattering) a more pronounced diffraction structure with sharp 
maxima and minima and a strong peak at zero angle as compared with the smoother angular 
dependence for other CD transitions. The increase of kinetic energy enhances asymmetry 
in the angular dependence of the /-averaged cross sections and decreases the role of the 
backward scattering (see Fig. 8). 



4.3 Hadronic atoms 

In order to illustrate the influence of the ns state energy shifts on the CD cross sections, 
we calculated the CD cross sections for (7rp) n + H collisions both with and without taking 
energy shifts into account. The effect is the most pronounced for the low- lying states and 
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Figure 8: The /-average dif- 
ferential cross sections for CD 
transition 6 — > 5 at different en- 
ergies. 
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is illustrated in Fig. 9. One can see that the maximal suppression due to the energy shift 
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Figure 9: The /-averaged CD cross sec- 
tions cr^j (thick lines) and a 1 ^ (thin 
lines) for (7rp) 4 + H collisions calculated 
with (solid lines) and without (dashed 
lines) taking the s-states energy shifts 
into account. The CD cross sections cr^ 3 
calculated both with (double-dot-dashed 
line) and without (dot-dashed line) tak- 
ing s-state energy shifts into account are 
also shown. 



of 4s state is about two times at very low energy both for 4 - 
while at E cm > 1 eV does not exceed 15% (for details see [14]). 
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Figure 10: The /-averaged CD cross sec- 
tions with An = 1,2,3 for the collisions 
of the pp atom (n = 8) with the hydro- 
gen atom. The dashed line shows the 
fit used in cascade calculations for the 
transitions with An = 1 and based on 
the mass-scaling of the results [6] for the 
muonic atom. The present results for 
the /-averaged elastic and Stark cross sec- 
tions are shown for comparison. 



The energy dependence of the CD cross sections for the collisions of the pp atom (n = 8) 
with the hydrogen atom obtained in the CC approach is shown in Fig. 10 for n = 8 and 
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the different values of An =1, 2 and 3. The special features of these cross sections are the 
following: the similar energy dependence but sharper than that of the elastic scattering and 
Stark transitions (see also in Fig. 10); the contribution of the transitions with An > 1 is 
comparable with the one for An =1 and is equal about 50%. The effect of the ns state shifts 
in the Z-averaged CD cross sections is negligible due to small statistical weight of the ns- 
state. In Fig. 10 we also compare our results with those obtained in the semiclassical model 
for the An = 1 transition. The satisfactory agreement is observed, but this agreement is 
quite occasional and takes no place for other n values. The distribution over the final states 
n' is completely different from the SC results [5] as it was illustrated in Fig. 10. 



5 Effect of closed channels 



The presented above results were obtained by the solution of the close- coupling equations 
including all the open channels. Although the channels with n > uq for uq < 7 and at 
energies less or ~ 1 eV are strongly closed, they can essentially change the open-channel 
wave functions at short ranges determining CD process. In the present studies we included 

i i i i i i i i i i i i 
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Figure 11: Dependence of the CD cross 
sections cr„ ^ no _i for (//~p) and (n~p) 
atoms at E cm = 1 eV on the closed chan- 
nels (n max is the maximal principal quan- 
tum number of the included channels). 



in the calculations all the open channels with n < n and step by step added closed 
channels with n' > n^, n' < n max to achieve the convergency of results. In Fig. 11 some 
of our preliminary results are shown for the Z-averaged CD cross sections at E cm = 1 eV. 
According to our investigation, the closed channel effect on CD cross sections depends 
crucially on the transition considered. As it is seen from Fig. 11 in case of (fi~p) atom 
the inclusion in the basis set of the closed channels results in a more pronounced effect for 
low-lying states. The investigations of convergency with the increase of number of closed 
channels are very time-consuming and are continuing now. 



6 Conclusion 

The unified treatment of the elastic scattering, Stark transitions and Coulomb deexcitation 
is presented within the quantum-mechanical close-coupling approach. The differential and 
integral cross sections for the above processes are calculated for the excited muonic, pionic 
and antiprotonic hydrogen atoms with n = 2 — 14 and relative energies relevant to the 
cascade calculations. The new results for CD process are obtained: anisotropy of the 
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angular distribution, substantial fraction of An > 1 transitions up to ~ 40% (for n > 4 at 
all energies under consideration), and a proper threshold behaviour of the CD cross-section. 
The calculated cross-sections are very important for the kinetics of the atomic cascade and 
give a more reliable theoretical input for the improved version of the cascade model [12] . 

We are grateful to L.Ponomarev, L. Simons, G. Korenman, T.Jensen and V.Markushin 
for fruitful discussions. This work was supported by Russian Foundation for Basic Research 
(grant No. 06-02-17156). 
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